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We investigate different types of neutrino hot dark matter 
with respect to structure formation and anisotropies in the 
cosmic microwave background radiation (CMBR). The possi- 
bility of neutrino hot dark matter produced through the decay 
of a heavier neutrino by the process ph — * vl + 4>, where is a 
scalar particle, is discussed in detail. This type of dark matter 
can possibly be distinguished observationally from the stan- 
dard neutrino dark matter by using new CMBR data from 
the upcoming satellite missions MAP and PLANCK. 
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The cold dark matter (CDM) model forms the ba- 
sis for all the currently viable models of the way struc- 
tures form in our universe. In its simplest form it con- 
sists of a critical (f2 = 1) density universe whose en- 
ergy density is almost solely (SIcdm — 0.95) in the form 
of non-relativistic non-interacting particles. In addition 
there is a small fraction of baryons (yis — 0.05). Struc- 
tures in this model grow from the original quantum fluc- 
tuations produced during the inflation era. The trou- 
ble with the simple CDM model is that, if the density 
fluctuation power spectrum is normalised to the CMBR 
quadrupole anisotropy measured by the COBE satellite, 
there are far too large fluctuations on small scales [Q. 
Any realistic model must therefore include something 
that dampens small scale fluctuations. Several possible 
solutions to this problem have been proposed, for exam- 
ple the inclusion of vacuum energy as a dominant compo- 
nent or tilting the initial power spectrum from inflation 
Perhaps the most promising model is the so-called 
cold+hot dark matter scenario (CHDM), where a compo- 
nent of hot dark matter (HDM) is included; particles that 
only become non-relativistic around the epoch of matter- 
radiation equality. Such light particles free stream and 
thereby smear out fluctuations below a certain length 
scale. This produces the desired result, namely that fluc- 
tuations on small scales should be lessened relative to 
the CDM model. A good fit to observational data can be 
obtained with Hhdm — 0.25 in a flat (Q = 1) universe. 

A light (eV) neutrino species seems by far the most 
realistic candidate for hot dark matter. If only one stan- 
dard neutrino species is available, its mass should be 
roughly 5-7 eV if Ohdm — 0.25 and ho ~ 0.5 (ho is 
the Hubble parameter in units of 100 km s _1 Mpc -1 ). 
However, numerical simulations have shown that a bet- 
ter fit can be obtained if there are two neutrino species 
available with half that mass (this model is referred to 



as 2J/CHDM); that is, a mass of roughly 2-3 eV g. This 
would demand a degenerate mass hierarchy or, alterna- 
tively, a four-component eV neutrino such as a light Dirac 
neutrino. Unfortunately, the sterile components of a light 
Dirac neutrino do not come into thermal equilibrium af- 
ter the QCD phase transition and therefore they only 
contribute very little to the cosmic energy density. A 
more plausible method of producing neutrino dark mat- 
ter is through the decay of a heavy neutrino to a final 
state containing light neutrinos, mimicking a degenerate 
mass hierarchy. 

There are several possible decay channels for a heavy 
neutrino to produce light neutrinos, for instance the ra- 
diative decay vjj — > ^L"f, the flavour violating weak pro- 
cess vjj — > v lvlvl or a decay of the type vh —* ^l4>, 
where <f> is a scalar or pseudo-scalar particle like the ma- 
joron |Q. Decays containing electromagnetically inter- 
acting particles are strongly constrained by observations 
so that in essence only decay to an "invisible" final state 
is allowed j5|. We shall look at the simplest example of 
such a decay, namely the vjj — » vl4> process. Specifi- 
cally we assume vh — v T and Vi, = v^. v e is assumed 
to be massless. A decay of this type occurs naturally in 
the so-called majoron models. In these models, neutri- 
nos acquire mass by interaction with a new Higgs field 
and the spontaneous breakdown of the global symme- 
try results in a Nambu-Goldstone boson called the ma- 
joron Q. In some of these models neutrinos have large 
Yukawa-type couplings to the majoron so that fast de- 
cays of the type vh — > vl4> are permitted. The coupling 
strength is very difficult to constrain experimentally, al- 
though a sufficiently large coupling would result in neu- 
trinoless double-beta decays with too short lifetimes [f|. 
For masses and lifetimes in the range we are interested 
in there are currently no direct experimental constraints. 

Neutrinos are in thermodynamic equilibrium with the 
cosmic plasma until the universe reaches a temperature 
of roughly 1 MeV, so that if the decay takes place prior 
to this freeze-out temperature, the additional density of 
light neutrinos will be washed out because of thermal 
equilibration. In the following we shall therefore always 
assume that the decay takes place subsequent to neu- 
trino decoupling, but prior to matter-radiation equality. 
In this case the number density of light neutrinos is twice 
that of a standard decoupled neutrino species, as in the 
2j/CHDM model. Note that if the decay takes place after 
decoupling, but at a temperature higher than T ~ 0.01 
MeV, the synthesis of light elements may be changed 
significantly, excluding some regions of mass-lifetime pa- 
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rameter space (6| . 

The purpose of this paper is to investigate observa- 
tional features of this "decaying neutrino cold+hot dark 
matter" (D^CHDM) scenario and compare it with the 
standard 2j/CHDM model. There are several differences 
between this decay scenario and the 2i/CHDM model. 
First of all, the distribution of light neutrinos will in 
general be non-thermal so that, even though the num- 
ber density of neutrinos in the two scenarios is equal, the 
energy density is different. Second, the amount of rela- 
tivistic energy density present will be different in the two 
cases. We will discuss the possibility that new CMBR 
experiments will allow us to distinguish between these 
different models of hot dark matter. In order to calcu- 
late CMBR and matter power spectra we have used the 
CMBFAST program developed by Seljak and Zaldarriaga 

0- 

The dynamics of the decay can be quite simply de- 
scribed with only two parameters, the mass, run, and 
lifetime, r, of the heavy neutrino. Using this we can 
write down the Boltzmann equations for the evolution of 
distribution functions 



df/dt-Hpdf/dp = C dec [f]. 
The decay terms are of the form M 



(1) 



||, a = 1/9 (m H /MeV) 2 (T H /s). This leaves us only one 
free parameter that controls the decay kinematics and 
thereby the shape of the particle distributions. The di- 
viding line between relativistic and non-relativistic decay 
is a = 1 so that the decay is relativistic if a < 1 and non- 



relativistic if 



> 



1. In Fig. 1 we show a realisation of 



this decay model for a = 1. 
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FIG. 1. The number density of different particle species 
as a function of temperature for a model with a = 1. As 
explained we assume vh = v T and ul = Up. v £ is assumed to 
be massless. no is the number density of a standard massless 
decoupled neutrino species. 



Cdec [f</>] = 
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where A(f H , f L , / ) = f„{\ - f L )(l + / ) - f L .U(l - f H ) 
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and TOq = m 

the lifetime of the heavy neutrino and g is the statistical 
weight of a given particle. We use gn = 9l = 2 and g^ — 
1, corresponding to <fi — <fi However, we could equally well 
have used g^ = 2 without changing the hnal results. Note 
that in all actual calculations we assume tol = to<£ = 
when solving the Boltzmann equation (this does not 
influence the results at all). When calculating transfer 
functions and CMBR power spectra we of course use the 
correct value for m^. The integration limits are 
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where i = L,<j>. We can vary the heavy neutrino mass 
without changing the final particle spectra if we keep con- 
stant the "relativity" parameter of the decay, defined as 
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FIG. 2. The effective temperature, pseudo-chemical poten- 
tial and amount of relativistic energy density in massless neu- 
trinos and scalar particles. To is the temperature of a stan- 
dard massless decoupled neutrino species. 
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Now, to simplify matters we parametrise the final light 
neutrino distribution in terms of an equilibrium Fermi- 
Dirac distribution with an effective temperature, T, and 
a pseudo-chemical potential, jit. In Fig. 2 we show these 
two parameters in units of the temperature of a stan- 
dard massless neutrino species. Further, we show the 
amount of relativistic energy density in massless neu- 
trinos (the energy density in the remaining third neu- 
trino species) and scalars in units of the energy den- 
sity of a standard neutrino species, iV„ ;re i. The energy 
density in light (HDM) neutrinos is not included here 
because this neutrino becomes non-relativistic around 
matter-radiation equality. Notice that for relativistic de- 
cays iVj, )re i is smaller than two because the energy den- 
sity in scalars is less than in a standard massless neu- 
trino species. For relativistic decays all these parame- 
ters approach constant values indicating that the decay 
proceeds in equilibrium, as expected. For non-relativistic 
decays, however, we see that the effective temperature in- 
creases while the pseudo-chemical potential drops. This 
reflects the fact that non-relativistic decays do not pro- 
ceed in thermodynamic equilibrium and the light neutri- 
nos are born with energies much above thermal. Also, 
the amount of relativistic energy density increases be- 
cause the scalar particles, <fi, are born with very high en- 
ergies. In fact, if the light neutrino mass approaches zero 
and the decay is non-relativistic this scenario approaches 
the rCDM model previously discussed in the literature 

To see how structure formation changes with different 
decay parameters we have calculated the matter power 
spectrum in terms of the quantity p2|, A 2 (k) = k , 
by using the CMBFAST package modified to include a 
non-equilibrium light neutrino species. We have then 
compared these spectra with the linear power spectrum 
data compiled by Peacock and Dodds |0|. Fig. 3 shows 
the power spectra for different models (note that all mod- 
els shown in Figs. 3 and 4 have been calculated with 
O = 1, Qb = 0.05, ho = 0.5 and spectral index n = 1). 
It is seen that for relativistic decays and equal f2„ our 
DfCHDM power spectrum closely resembles that of the 
2j/CHDM scenario. This is not too surprising because 
the light neutrino distribution is close to thermal in this 
case. For non-relativistic decay, the difference is pro- 
nounced. The power spectrum breaks away from that 
of the CDM model at lower fc-values if a 3> 1 than if 
a < 1 because of the larger free-streaming length of the 
light neutrinos. Even at large k the power spectrum has 
less power if a is large. This is because the amount of 
relativistic energy density is significantly larger if the de- 
cay is non-relativistic. The figure also shows that varying 
fi„ influences the power spectrum differently than chang- 
ing a, because varying Q v changes the amount of CDM 
and the free-streaming length of light neutrinos, whereas 
varying a changes the amount of relativistic energy den- 



sity and the free-streaming, but not the amount of CDM. 
We have also shown the limiting case of large a and low 
where our model closely resembles the rCDM model. 
As already shown in Ref. [ flO| this model provides a very 
good fit to the linear power spectrum. 
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FIG. 3. Matter power spectra with HDM stemming from 
neutrino decay. The two sets of lines correspond to a = 0.1 
and 10, the upper curves being a = 0.1 and the lower a = 10. 
For comparison we also show the 2^CHDM power spectrum, 
as well as the rCDM-like power spectrum with S7„ = 0.01, 
a — 70 and ho = 0.45. The diamonds represent the compiled 
linear data of Peacock and Dodds |l2| . 

Next, we also look at the power spectrum of CMBR 
anisotropies produced by these scenarios. It has been 
known for some time that the CHDM scenarios produce 
larger CMBR fluctuations than the CDM model and that 
the peaks in the power spectrum are also shifted Q . We 
have calculated the power spectra for our series of mod- 
els in terms of C; coefficients, C\ = (|a; m | 2 ), where the 
ai m coefficients are defined in terms of the temperature 
fluctuations as T(0, <f) = J2i m a lmYlm(0, 4>)- 

Fig. 4 shows power spectra for several different models, 
illustrating the differences between the CDM, 2z/CHDM 
and D^CHDM models. Exactly as with the matter trans- 
fer function we see that for relativistic decays, the CMBR 
power spectrum for the D^CHDM model is very similar 
to that for the 2^CHDM scenario, indicating that the de- 
cay proceeds roughly in equilibrium. For non-relativistic 
decays, however, the power spectrum for the D^CHDM 
model is markedly different, especially the amplitude of 
the first Doppler-peak is much larger. The reason is that 
increasing a corresponds to increasing the relativistic en- 
ergy density, meaning that recombination occurs much 
closer to the radiation dominated epoch than in the stan- 
dard model. As explained for instance in Refs. [ ]l0| , |l4| |, 
the reason is that in the matter dominated epoch pho- 
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tons travel through a constant gravitational potential af- 
ter they are emitted from the last scattering surface. In 
the radiation dominated epoch this is not the case and if 
the universe is still close to being radiation dominated at 
recombination the non-constant gravitational potential 
increases the CMBR anisotropy, a feature known as the 
early ISW effect pl| . Changing the HDM content by rea- 
sonable amounts in the different models only influences 
the CMBR power spectra slightly O. 
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FIG. 4. CMBR power spectra for CDM, 2^CHDM and two 
different D^CHDM models. All the power spectra have been 
normalised to the quadrupole coefficient 6C2. 

In conclusion, we have seen that a natural way of pro- 
ducing a 2^CHDM like model is by decay of heavy neutri- 
nos. If the decay proceeds relativistically this Dj/CHDM 
scenario closely resembles the standard 2^CHDM sce- 
nario because particle distributions remain roughly ther- 
mal. For non-relativistic decays there can be large dif- 
ferences, both in the matter power spectra and in the 
CMBR fluctuation spectra, compared with the 2z/CHDM 
model. We have also shown that limiting cases for our 
series of models is, in one end, the rCDM model [[Toll , 
where —> and -/V„. rc i is large (non-relativistic neu- 
trino decays with very low light neutrino mass), and in 
the other extreme the 2^CHDM model (relativistic de- 
cay with light neutrino mass equivalent to that in the 
2^CHDM model) |3]. A decaying r neutrino can there- 
fore change structure formation on both small and inter- 
mediate scales. Changing the time of matter-radiation 
equality changes the amount of structure produced on 
horison-size scales at this time, whereas changing the 
amount of neutrino hot dark matter changes the amount 
of small scale structure. There is some real hope that 
within the next few years it will become possible to dis- 
tinguish between these different models. Both because 
new and improved galaxy surveys like the Sloan Digital 



Sky Survey p6| are forthcoming and should greatly im- 
prove our knowledge of the matter power spectrum and 
because the next few years will hopefully see a leap in 
our understanding of the cosmic background radiation 
because new missions are scheduled to fly shortly, mea- 
suring the angular power spectrum to I — 1000 and be- 
yond to great precision [jljj . Using this new data it should 
be possible to constrain structure formation models like 
the Dj/CHDM model described here very strongly. 
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